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Abstract: We discuss the computation of correlation functions in holographic RG flows. 
The method utilizes a recently developed Hamiltonian version of holographic renormaliza- 
tion and it is more efficient than previous methods. A significant simplification concerns 
the treatment of infinities: instead of performing a general analysis of counterterms, we 
develop a method where only the contribution of counterterms to any given correlators 
needs to be computed. For instance, the computation of renormalized 2-point functions 
requires only an analysis at the linearized level. We illustrate the method by discussing flat 
and AdS-sliced domain walls. In particular, we discuss correlation functions of the Janus 
solution, a recently discovered non-super symmetric but stable AdS-sliced domain wall. 
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1. Introduction 

One of the successes of the gravity/gauge theory correspondence is the ability to perform 
detailed computations of QFT correlation functions by doing a classical gravitational com- 
putation. The holographic prescription put forward in jl], || identifies the boundary values 
of bulk fields with sources of gauge invariant operators of the boundary theory and the 
bulk partition function (which is a functional of the boundary data) with the generating 
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functional of correlation functions of gauge invariant operators. In particular, the super- 
gravity on-shell action is identified with the generating functional of connected graphs of 
the (strongly coupled) boundary QFT at large N. 

The correlation functions computed using this prescription exhibit divergences due to 
the infinite volume of spacetime and need to be renormalized. A systematic method for 
doing this was developed in [j^, ||, |5| following earlier work in ||, see Q for a review. The 
method is complete and can be used to compute renormalized correlation functions in any 
holographic RG-flow described by a (locally) asymptotically AdS (AAdS) spacetime. Some 
of the steps, however, are rather tedious which make the method difficult to apply. In this 
paper we will present a significantly simpler version of the method. 

Let us first recall the computation of 2-point functions as developed in the early period 
of the AdS/CFT correspondence ||. We will call this method the "old approach". To 
regulate the theory one imposes a cutoff at large radius and solves the linearized fluctuation 
equation with Dirichlet boundary condition at the cutoff. The second variation of the cutoff 
on-shell action is then computed in momentum space yielding an expression containing 
singular powers of the cutoff times integer powers of p plus non-singular terms in the 
cutoff which are non-analytic in p. The 2-point function is defined as the leading non- 
analytic term. The polynomial terms in p which are dropped are contact (^-function terms 
in the position space correlator, which are scheme dependent in field theory and largely 
unphysical. The non-analytic term has an absorptive part in p which correctly gives the 
2-point function for separated points in x-space. This method is quite efficient, but it 
is not fully satisfactory since it is not correct in general to simply drop divergent terms 
in correlation functions. The subtractions should be consistent with each other and they 
should respect all (non-anomalous) symmetries. 

In holographic renormalization one replaces the above computation by a two step 
procedure. In the first step renormalization is performed. This is done via the so-called 
near-boundary analysis. This step associates to a given bulk action a set of universal local 
boundary counterterms that render the on-shell action finite on an arbitrary solution of the 
bulk equations of motion. Furthermore, to each bulk field we associate a pair of conjugate 
variables: the "source" and the exact 1-point function in the presence of sources (or "vev"). 
These can be read off from the asymptotic expansion of the bulk fields. 

n-point functions can now be obtained by functionally differentiating the 1-point func- 
tion w.r.t. the sources. The near-boundary analysis leaves undetermined the vev part of 
the solution. To obtain the vev part in terms of the source part we need an exact (as 
opposed to asymptotic) solution of the bulk field equations with arbitrary Dirichlet data. 
Such computation is (presently) not possible in all generality but one can solve the bulk 
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equations perturbatively around a given solution. In particular, to obtain 2-point func- 
tions it is sufficient to obtain a solution of the linearized field equations, thus connecting 
with the method described earlier. Since all subtractions are made by means of covariant 
counterterms mutual consistency is guaranteed. Extra bonus of the method is that all 
Ward identities and anomalies are manifest ab initio and 1-point functions (with sources 
equal to zero) are automatically computed. Furthermore, the procedure demonstrates that 
holography exhibits general field theoretic features such as the fact that the cancellation 
of UV divergences does not depend on the IR physics. The main drawback of the method 
is that the near-boundary analysis is somewhat tedious to carry out. 

It is conceptually satisfying that the procedure of renormalization can be carried out 
in full generality without reference to a particular solution, i.e. that the counterterms 
are associated to a bulk action, not a particular solution. On the other hand, however, 
interesting solutions corresponding to different RG-flows usually involve different sets of 
fields and different actions. To compute correlation functions in different RG-flows one 
would thus need to first complete the near-boundary analysis for each different action. 
Furthermore, some counterterms may have vanishing contribution when we consider a 
specific solution and may contribute to only a few correlation functions. Thus, if we are 
only interested in the computation of a small set of correlation functions it would be more 
efficient if we had a method that computes only the contribution of the counterterms to 
these correlators (rather than first computing all counterterms and then specializing to the 
specific case). The main result of this work is the development of such a method. 

In ||, extending previous work |u], [HJ 12 1 (see also [13|), we developed a Hamiltonian 
method to obtain renormalized correlation functions of the dual quantum field theory from 
classical AdS gravity. The method is based on a Hamiltonian treatment of gravity using the 
AdS radial coordinate as the 'time' coordinate. More concretely, the canonical momenta 
conjugate to the bulk fields evaluated on a regulating hyper surface are identified with the 
(regulated) one point functions. These are functionals of the bulk fields on the same hy- 
persurface and the bulk field equations become first order functional differential equations 
for the momenta. To obtain covariant counterterms one expresses the radial derivative in 
terms of functional derivatives w.r.t. the bulk fields. The Dirichlet boundary conditions 
on the bulk fields then imply that the radial derivative is asymptotically identified with 
the total dilatation operator of the dual field theory. This is indeed consistent with the 
interpretation of the bulk radial coordinate as the energy scale of the dual field theory. The 
momenta are then expanded in covariant eigenfunctions of the dilatation operator, which 
is the analogue of the asymptotic expansion of the bulk fields in the standard approach. 
However, since the expansion is by construction covariant, all counterterms, including those 
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related to the conformal anomaly, are directly computed without the need to invert the 
asymptotic expansions. The method, as the standard method of holographic renormaliza- 
tion, makes all (anomalous and non-anomalous) Ward identities manifest and in fact we 
showed that the two approaches are completely equivalent. Nevertheless, the new method 
leads to a much more efficient algorithm for the computation of correlation functions from 
holography as we will demonstrate in this paper by various examples. The main advantage 
is that the divergences of the holographic correlators can be removed consistently without 
the need for a general near boundary analysis as was done in previous methods. 

The paper is organized as follows. In Section 1 we recall some of the highlights of the 
method developed in Q to the extent that is necessary for the subsequent discussion. In 
Section 2, we consider Poincare domain wall solutions of the bulk supergravity theory with a 
single active scalar field turned on. We discuss how to distinguish between explicit breaking 
of the conformal symmetry by an operator deformation and spontaneous breaking by vevs 
by computing the exact one-point functions for these backgrounds. Two point functions 
are then calculated and we show explicitly how they can be renormalized without the need 
for a general near boundary analysis. Finally, in Section 4 we study the recently found 
Janus solution [14], which is a particular non-supersymmetric but stable AdS domain wall 
1 15]. We calculate the vevs of the background as well as some two-point functions and 
we show that the Ward identities associated with the symmetries of the background are 
satisfied. 

2. Hamiltonian Holographic Renormalization 

We start by considering the truncated supergravity action describing gravity in (d+1)- 
dimensions coupled to a single scalar field: 

%, <&] = f d d+1 x^g \~^R + + V(Q)] - i / d d x^2K. (2.1) 

Jm y IK I J IK J dM 

where dM is the conformal boundary of the (locally) asymptotically AdS spacetime M. 
By choosing a suitable gauge the bulk metric can be taken to be of the form 

ds 2 = g lxv dx tJl dx v = dr 2 + 7y (r, x)dx l dx j , (2.2) 

where i = 1, . . . d. Using the Gauss-Codazzi relations the bulk equations of motion are 
reduced to Si 



K 2 - KaK lJ = R + k z 



K) + KK) = R) - k 



$2 _ yj <Q.$g.$ _ 2V($) , 

Vi-RTj - VjK = K 2 <S>dj$, (2.3) 

d^dj® + - — : 7 V(<f>)5) 
d — 1 J 
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Here K^j = ^jij is the extrinsic curvature and K l - stands for ^(7 K^j). Additionally, we 
have the equation of motion for the scalar field 

<t> + K<f> + \3<S>-V'{$) = (2.4) 

and the on-shell action takes the form || 

Son-shell = d d X^(K-X), (2.5) 



where A satisfies the equation 



Ik 2 

\ + K\ + j-jV{$) = 0. (2.6) 



The Hamiltonian formalism allows one to write the canonical momenta on a regulating 
surface S r , diffeomorphic to the boundary dM at r — > 00, as functional derivatives of the 
on-shell action on S r w.r.t. the corresponding bulk field. This shows that the canonical 
momenta (on S r ) are covariant functionals of the induced fields on S r , which in turn leads 
to an expression for the radial derivative in terms of functional derivatives: 

d r = J dfiaaKijh,*]-^- + J d d x^,^. (2.7) 

The asymptotic behavior of the bulk fields then implies that as r — > 00 

8 r ~ I d d x2 li ~ + f d d x(A - d)<S>^ = 5 D , (2.8) 

J hij J 0$ 

which is precisely the dilatation operator. This observation is at the center of our approach. 

The next step is to expand the canonical momenta and the on-shell action in eigen- 
functions of the dilatation operator as 

Ki[ 7 , d>] = K {0) ) + K m ) + ■■■ + K (d) ) + K (d) ) loge" 2 ' + • • • , 

A[ 7 , $] = A ( o) + A {2) + • • • + \ {d) + \ d) log e" 2r + • • • , (2.9) 

?r [ 7; = = ^ 7T W + 7T (A) + 7T (A) log e" 2r + • • • , 



d-A<s<A 



where 



S D K (n)j = ~ 


-n K (njj> n<d, SdK^j = - 


-dK (d) ), 




5 D K (d)) = -dK {d) ) - 






) = -nA (n) , n < d, S D \d) = 


—d\d): 




^D^(d) = -d\d) 


~ ^\d): 


^D^(s) = -S7T(s): 


, d — A < s < A, 5dtt(a) = " 








- 2tt(A) 



(2.10) 
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In particular, the coefficients of the logarithmic terms are related to the conformal anomaly 1 , 
while the terms which transform inhomogeneously under the dilatation operator are in gen- 
eral non-local and correspond to the exact one-point functions of the dual operator. More 
precisely, we have 

(Tij)ren = J ( K (d)ij ~ K {d)lij) , (O)ren = ~p 7T(A)- ( 2 -ll) 

K vT 
At the same time, the renormalized on-shell action is given by 

Sren = ~ d d X^j{K (d) - \ d) ). (2.12) 

The second equation in fl2.3|), then leads to the Ward identity 

Vi(lj} r en = -(0} ren d^, (2.13) 

whereas the following relation, which was derived in ||, 

(1 + 5 D )K (d) + k\A - d)vr (A) $ = {d+ 5 D )\ d) (2.14) 
gives the (anomalous) Ward identity 

(Tt) Tcn = (A - d)<Z>(0) rcn + A. (2.15) 
Here A = —-^(Kr d \ — \(d)) is the conformal anomaly defined by 

S D Sren = - dtxyfiA (2.16) 

As a consistency check, notice that the above transformation rules ( |2.10| ) imply that the 
conformal anomaly is conformally invariant as required: 

S 2 D S rcn = 0. (2.17) 

A general algorithm for the recursive evaluation of the coefficients in the covariant 
expansions was given in Q. It leads to general expressions for covariant counterterms in 
arbitrary dimension for pure gravity, but in general the details depend on the matter that 
couples to gravity and so we will skip the details except when necessary to illustrate the 
examples we give below. A final comment is due regarding the renormalization scheme 
dependence of the correlators we have computed. One could add extra finite covariant 
boundary terms in the action ( [2.1| ) without altering the equations of motion. The corre- 
sponding terms in the covariant expansions would then be shifted accordingly and so would 
be the one-point functions. This is all in complete agreement with our expectations from 

the field theory point of view. 

1 Recall that the gravitational part of the conformal anomaly is non-zero only for even dimension d but 
matter conformal anomalies jl6) are generically non-zero for all d. Correspondingly, the logarithmic terms 
for pure gravity are non-zero only for even boundary dimension d, but they are generically non-zero even 
for odd d when additional fields are turned on pj. 



-6- 



3. Poincare Domain Walls 



In this section we will consider linear fluctuations around Poincare domain wall solutions 
of the equations of motion. These take the generic form 



ds% = dr 2 + e 2A{r) dx i dx i , $ = <j> B (r 



(3.1) 



Inserting this ansatz into the equations of motion we find that A(r) and <^_B(r) satisfy 



d(d - 1) 



^ - 2V(0 B ) = 0, 



A + dA 2 + 



2k 1 



-V(<t> B ) = 0, 



d-1 

B + dA<j) B - v'{<p B ) = 0. 



Moreover, one has 



which immediately implies 



\ B + dAX B + 



2k 2 
d-1 



V^b) = 0, 



X B = A + £(x)e 



-dA 



-K B + 



(3.2) 

(3.3) 
(3.4) 

(3.5) 
(3.6) 



d 

where ^ is an arbitrary integration function of the transverse coordinates. Therefore, 

Sfn-sheii = — J-?: \ d d x^K B + ^ d d x£(x). (3.7) 

The last term corresponds to finite local counterterms. As expected, there is an ambiguity 
on the on-shell value of the action corresponding to the renormalization scheme dependence 
of the dual field theory. 

It is well known that the second order flat domain wall equations are solved by any 
solution of the first order flow equations 

.2 



A 



K 



-W{<t> 



B 



d-l 



provided the potential can be written in the form 



V(<f> 



B) 



W '2 _ ^L W 2 
d-l 



(3.8) 



(3.9) 



The motivation for considering theories with such potentials stems from the fact that they 
guarantee gravitational stability of the AdS critical point and of associated domain-wall 
spacetimes, provided the AdS critical point is also a critical point of W [17, 18, 15]. This 
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means, by the AdS/CFT duality, that the dual theory is unitary. Notice that in principle it 
is always possible to write the potential in the form ( |3.9D if one views ( |3.9[ ) as a differential 
equation for W(<Ab). The resulting W however may not have the original AdS spacetime 
as a critical point. Furthermore, in practice it is considerably difficult to solve ( |3.9| ). As a 
first step in this direction we observe that (3.9) can be transformed into Abel's equation Q. 
In general, the set of solutions of the second order equations may include solutions which 
cannot be obtained from the first order flow equations. In this section we will restrict 
attention to solutions which can be derived from the flow equations. For this class of 
solutions we have 

Son-shell = [ d d x^W(<f> B ) + -^ [ d d x£(x). (3.10) 
3.1 Deformations vs VEVs 

The domain wall solutions we have described correspond via the AdS/CFT duality to defor- 
mations of the boundary CFT by relevant operators, in which case the conformal invariance 
of the boundary theory is explicitly broken, or to a vacuum expectation value of a scalar 
operator which spontaneously breaks the conformal symmetry. Marginal deformations are 
also similarly described. 

Locally asymptotically AdS metrics satisfy the asymptotic condition 7ij(r, x) ~ e 2r ^/ij^(x) 
as r —* oo, which requires A(r) ~ r. Moreover, a scalar field dual to an operator of dimen- 
sion A behaves asymptotically as $>(r, x) ~ e~^ d ~^ r (j)^(x). These asymptotic conditions, 
together with the equations of motion, require that the scalar potential takes the form 

V(*) = -^|^ + W + (3.11) 

where the mass is related to the dimension A of the dual operator by m 2 = A(A — d). 
Solving the equation of motion for the scalar field with such a potential leads to a generic 
solution of the form 

*(r,x) = e-( d ~ A > [0 (o) (x) + ...]+ e~ Ar [<P {2A . d) (x) + ...]. (3.12) 

We will consider operators for which d — A > A, or A > d/2. Moreover, we are interested 
in relevant or marginal operators and so A < d. In total then d/2 < A < d. In this 
range the first term in the solution for <I> is dominant asymptotically and corresponds to 
the source, while the second term is related to the one-point function of the dual operator. 
In the special case A = d/2, which saturates the BF bound, the scalar field takes the form 2 

$(r, x) = e~ dr / 2 \-2r (0 (o) (*) + ...)+ (o) (x) + . . .1 . (3.13) 



2 Note that the radial coordinate we use here is related to the Fefferman- Graham radial coordinate by 
p = e~ 2r . The factor — 2r = logp is chosen to match with the corresponding formulae in W. H. 
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Again, the first term is the source for the dual operator, while the second is related to its 
expectation value. 

Depending on the form of the 'superpotential' W{4>b) the domain wall solution can 
describe either a deformation of the dual CFT or a phase with spontaneously broken 



conformal symmetry as we now explain. A similar analysis can be found in [12|. Assuming 
the potential has a critical point at 4>b = 0, equation ( |3.9| ) together with the flow equations 
Q3.8D and the requirement that 4>b = is also a critical point of W, imply that W has an 
expansion around (j>B = of two possible forms: 

W + {4> B ) = - \(d - A)0| + . . . (3.14) 



W_(0b) = -^-^A^ + .... (3.15) 



Which of these cases is realized is purely a property of the background solution and we 
need to examine each case separately. Moreover, the extremal values d/2 and d of the 
scaling dimension A require special attention. We will now compute the vev for the stress 
energy tensor and dual scalar operator for all cases. 

As was shown in || , the part of the divergent part of the on-shell action involving only 
the scalar field is a function U($>), satisfying the equation for the 'superpotential' W (p^Sj ), 
and having an expansion around $ = 

U($) = -^i_I(d-A)<I> 2 + ... (3.16) 



K 



Let us consider first the case W+ is realized in the background. In this case we can 
choose a scheme where the counterterm action is W+. To see this notice that any two 
solutions of ( |3.9| ) with identical expansions around <3? = up to order $ 2 can only differ at 
order ~ e~ dr . This is easily proved by looking for the most general power series 3 

solution of ( |3.9| ) with this particular form up to quadratic order in <£. One finds that all 
terms are uniquely determined up to order d/{d — A) where the recursion relations break 
down. This is precisely where an arbitrary integration constant appears and the two solu- 
tions could be potentially different. However, this is irrelevant for the purpose of removing 
the divergences of the on-shell action and so we can choose the renormalization scheme 
J7(3>) = W+(<I>), and set the integration function £ to zero. This choice of counterterms 
corresponds to a supersymmetric renormalization scheme since it ensures S^ n = Q. It 
follows that the background vevs of both the operator dual to and the stress tensor 
vanish identically and so this background describes a deformation of the boundary CFT 
by a relevant operator. 



3 As usual one must include a logarithmic term at order d/(d — A) in the general case. 
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Let us now consider the case W- is realized in the background. In this case we cannot 
choose W- as the counterterm since it differs from U(&) at the quadratic order. In this 
case, however, 4>b ~ e~ Ar and so the on-shell action evaluated on the background contains 
only the volume divergence since, by the hypothesis, 2A > d. Hence, again, setting £ = 
corresponds to a supersymmetric renormalization scheme with 5^ n = 0. Accordingly, 
the background expectation value of the stress tensor vanishes, but not that of the scalar 
operator. In this case W'_{4>b) — U'(4>b) = — (2A — d)(f>B + • ••, and hence, 

(0)g n = (d-2A)<f> B , (3.17) 

which spontaneously breaks the conformal symmetry of the dual CFT. Here we used the 
fact that the regularized 1-point function (=canonical momentum) is related to the super- 



potential via the first order equation (3^). The renormalized 1-point function is obtained 
by subtracting the contribution U' of the counterterm. 

It remains to examine the two extremal cases A = d/2 and A = d. When A = d/2 
there is no distinction between W+ and W— as they are equal and 4>b ~ e~ dr / 2 , i.e. it 



behaves asymptotically as the vev term in (3.13). However, the on-shell action function 
[/($) includes divergences coming from the source term in ( 3.13; ) and therefore it cannot 
be identified with VF( ( I ) ). It is straightforward to find the covariant counterterms for this 
case with the method of ||] , but the singularity structure for the terms involving the scalar 
field are not the standard ones, as we now explain. This can be traced back to the fact 
that the source term for the scalar contains a logarithm, in contrast to the generic case. 
A simple calculation using the asymptotic form of the solution shows that the canonical 
momentum of the scalar field takes the form 



vr = $= (- -- )$ + ... (3.18) 
r 2 



Hence, the on-shell action is 



Son-shell = ^_ shell + jf d d X^\ ( I - ^ $ 2 + . . . (3.19) 

where S^n-sheii is the on-shell action for pure gravity. We have therefore shown that 



d — 1 1/1 d \ , 
tf(*) = _ + -( --)** + ... (3.20) 



which is the divergent part of the on-shell action and must be removed (of course there is an 
other part coming from pure gravity). The same counterterms (for d = 4) were derived in 
[J]. ||. We now see that the scalar operator gets a vev since W'(^>b) — U'((/)b) = — \4>B + • • 
i.e. 

(0}g n = 2ct>B- (3.21) 
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This also agrees with the results in j|, ||. Again the difference between W(<Pb) and U(4>b) 
is subleading and the stress tensor gets no vev since S ren — 0. 
Finally we consider the case A = d, for which 

$(r,x) = [<P (0) (x) + ...]+ e~ dr [<P {d) {x) + ...]. (3.22) 

The equations of motion require V'(<3?) = and so the potential is just the cosmological 
constant V(&) = — ^4-^ . It follows that the on-shell function U($>) is also a constant, 
i.e. the first term of W±. In this case, however, the general solution to ( |3.9D can be easily 
obtained. There are two distinct solutions (cf. eq. (2.10) in |l5| ), 

W+ = -^±, W_ = -^co8h^^fo-^ (3-23) 

in agreement with our general analysis. Notice that in this case the supergravity action 
and hence the second order field equations are invariant under constant shifts of the scalar 
field. Such a constant corresponds to the source term of the solution. One may use this 
symmetry to set cfto to zero in W-. Then, exactly as for the case d/2 < A < d, if W+ is 
realized in the background, then neither the scalar operator nor the stress tensor acquire a 
vev, and if W_ is realized, the scalar operator gets a vev —dcfts, while the vev of the stress 
energy tensor vanishes. 

So finally we can summarize all possibilities for flat domain wall backgrounds in the 
following table: 



A 


W 


(O)fen 


(Ti) B 

\^ j / ren 


d/2 < A < d 


+ 












(d - 2A)4> B 





d/2 


± 


2(f>B 






3.2 Linearized Equations 

Now let us consider fluctuations around the backgrounds we have described so far. We will 
only keep terms up to linear order in fluctuations, which suffices for the calculation of the 
two point functions. The metric fluctuations take the form 

lij = lij (?) + hij (r, x) = e 2A ^5ij + (r, x) , (3.24) 

and the scalar field is 

$ = <Mr) + <Kr,z). (3-25) 
The extrinsic curvature then becomes 

K} = A5] + ^S}, (3.26) 
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where 5j = Jgh^j. Sj can be decomposed into irreducible components as 

,/ / 1 ;)'/) \ ;)';)• 



5' = e' + <f £j - + V + — y- /s - J , + _^ s, (3.27) 

where 9je*- = = die 1 = and indices are raised with the inverse background metric 
e -2Agij Conversely, each of the irreducible components can be expressed uniquely in 



terms of Sj as 



e} = IT fe jSi , e i =7r i —S l , f = ir k S l , 5 = 4^, (3.28) 
where we have introduced the projection operators 



ny, = ( *w, + T a T*i - 7— ) > ( 3 - 29 ) 



2 \ j ^ d - i 

and 

^ = 65-— 2. (3.30) 

With this nomenclature we can now go on and derive the equations of motion for the linear 
fluctuations. The result is 4 

(d? + dAd r + e" 2A n) e) = 0, (3.31) 
(d 2 T + [di + 2W&1 log VK]d r + e~ 2A dj u = 0, (3.32) 

/ = -2k 2 <£b0, (3.33) 



S 



1 



(d- l)i4 L 



e □/ + 2/t ( (j>B(t> — V{(j)B)(j)} (3.34) 



where 



and we have used the diffeomorphism invariance in the transverse space to set e% = 0. The 
last two equations give immediately the momenta dual to / and S and hence the corre- 
sponding one-point functions with linear sources. Moreover, since the canonical momenta 
are functionals of the bulk fields || , to linear order in the fluctuations we must have 

e) = E(A, <f> B )e), co = Q(A, <j> B )u. (3.36) 

The first two equations then become first order equations for E and f2: 

E + E 2 + dAE - e" 2 V = 0, 

n + n 2 + [dA + 2W&1 log w}n- e~ 2A P 2 = o, (3.37) 



4 Note U B = e~ 2A U = e- 2A & t3 didj 
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where we have performed a Fourier transform in the transverse space. Given the solutions 
for E and O we can immediately write down all momenta, namely 

e) = Ee), (3.38) 
/ = -2K 2 <j) B cf>, (3.39) 

= (W + n)0 + ^L^n/, (3.40) 

To completely determine the one-point functions with linear sources we first need to obtain 
exact solutions for E and O and secondly, to determine the covariant counterterms for the 
momenta, but only to linear order in the fluctuations. Since e~ 2A and W(4>b) are already 
covariant functions of the background fields, it suffices to find covariant expansions for 
E and £1 in the background fields. These can be organized according to the dilatation 
operator for the background 

5 D = d A + (A-d)<f ) Bd< t , B . (3.42) 

More generally, the radial derivative is expanded in functional derivatives w.r.t. the back- 
ground fields as 

K 2 

8 r = Ad A + fl fy B = -j—W(<p B )dA + W'{<p B )d <t>B ~ 8 D + . . . (3.43) 
Inserting the following expansions 5 for E and Q in the first order equations (|3.37| ), 

E = E (l) + ■■■ + E (d) log(e- 2r ) + E (d) + • • • , 

n = n (0) + ■■■ + n {2A _ d) io g ( e - 2r ) + n (2A _ d) + ■■■ , (3.44) 

one determines all covariant counterterms which render all momenta finite to linear or- 
der in the sources. This procedure is substantially simpler than the general holographic 
renormalization required to determine the full non-linear counterterms and is a significant 
improvement over previous methods. 

A final simplification can be made for the case of backgrounds corresponding to defor- 
mations of the dual CFT. As we saw in the previous section, the 'superpotential' W of the 
background can be included in the counterterm action, corresponding to a supersymmetric 
renormalization scheme. After this counterterm is added to the on-shell action (we still 

5 These expansions are strictly correct for d/2 < A < d, but we will deal with the special case A = d/2 
in the examples below. 
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have to determine the counterterms for E and f2), the momenta take the simpler form 

Ee), (3.45) 



S = ~- 



/ = o, 
i w 

2k 2 W 

w 



-2 A 



w 



-□ 



(3.46) 
(3.47) 

(3.48) 



3.3 Examples 

We will treat the two examples that have been the main testing ground for holographic 
computation of correlation functions, namely the GPPZ flow [19] and the Coulomb branch 
flow [20, pi] . The computation of certain two-point functions for the CB flow was first 



discussed in j2fl] and for the GPPZ flow in [22]. Two-point functions for both flows were 
systematically studied in ]^3|, |4j], see also ]^4|, 25, ^] for earlier work. Since the results are 
known, the emphasis here will be in method rather than the correlators themselves. 
GPPZ Flow 

The GPPZ flow describes a deformation by a supersymmeric mass term of N = 4 
SYM. The bulk theory is that of a scalar field dual to an operator of dimension A = 3 
coupled to gravity in five dimensions. The background 'superpotential' is 



W(cPb) 



3 
2^2 



1 + cosh 




(3.49) 



which is of the form W+ and corresponds to a deformation of the boundary CFT by a 
relevant operator. The background solution takes the form 



>B 



~J 3 2 l0g 



u 



l-VT 

It is also useful to note the relations 



u 



„2A 



1 



1 



-2r 



U 



W 



_3_1 

9 ' 



w' 



u 



W" 



2k 2 



Changing variable from r to u in ( 3.37| ) we obtain 



2(1- u)E'(u) + E 2 + -E 
u 



2(1 - u)Q' '(it) + Q 2 



2 n 



u 

l-u 



-W + l. 

0, 
= 0. 



II 



(3.50) 
(3.51) 

(3.52) 
(3.53) 



The solutions which are regular at u = are 



E{u) 



-/(l-u) 



F 1 



2,1+ 2° ! 3; u 



F 



2i 2 ' ' 



(3.54) 
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and 

1 TP p-g 3+a . o. \ 

Q(u) = -p 2 (l-u) 2 \l ' ' (3.55) 



where a = \J\ — p 2 . 

Next we need to find covariant counterterms for E and 0. Inserting 



d r = 5 D + —<j) B ( d A - ^Bd^g ) H (3.56) 



(3.57) 



and the expansions ( [3.44 ) in ( |3.37| ) one very easily determines 

771 P 2 -2A i P 2 -2A ( P 2 -2A i ^ 2 1 2 \ i -2r , 77. 

E = Y e + T e IT 6 + y^ b ) ge + ■%) + •■•. 

2 

= -ye- 2A loge- 2r + (2) + --- . (3.58) 

Expanding the exact solution in 1 — u and removing the covariant terms we have just 
determined allows for the evaluation of Eu\ and f2( 2 ), which are precisely the terms required 
to calculate the renormalized one-point functions. Putting everything together we find the 
following two-point functions: 

(o(p)o(- P )) = -~p 2 j, (Ti(p)o(-p)) = y^p'j, 



(Ti(p)Ti(-p)) = -_p 2 (J + 1), p> Pi (Ti) = 



3 

n 2 

(T l3 (p)T kl (-p)) TT = ^Il ljkl [±p 2 (p 2 + A)K + ^) 



where 



J = 2^(1) - ^ + \^V 2 ) - - 
^ = ^(l)+V(3)-^(2 + |)-V(2-|) 
Coulomb Branch Flow 

The Coulomb branch flow is a solution of five dimensional AdS gravity coupled to 
a scalar field of mass m 2 = —4, which therefore saturates the BF bound. The solution 
describes the case where an operator of dimension 2 gets a vev. The superpotential is 



W(<j> B ) = -\ 



(3.59) 



The solution can be parametrized by v = e^ K ^ B as 

r = 2v 2/3 (l -v), e~ 2A = v~ 2/3 (l-v), W = --^v^ l ^(v + 2). (3.60) 
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The boundary is located at v = 1. In terms of v the first order equations ( 3.37 ) become 



2(1 _ v)E'(v) + v- 2 / 3 E 2 + i (l + ^ E - pV 4 / 3 (l 



v)=0, 



and 



2(i-v)n'(v)+ v - 2 / 3 n 2 + Ui + 



12 

v + 2 



Q-p 2 v- 4/3 (l-v) = 0. 



The solution for E which is regular at v = is 



E(v) = 2a(l - v)v 



-1/3 



1 + 



av F{a + 1, a + 1; 2a + 3; u) 
2(a + l) F(a, a; 2a + 2; v) 



(3.61) 



(3.62) 



(3.63) 



where a = — | + ^y/l + p 2 . We will not give here explicitly the exact solution for Q since 
it is rather complicated. To obtain such a closed form solution one must transform the 
above equation for £1 into a soluble form and then obtain Q implicitly through the solution 
of the transformed equation. After obtaining covariant counterterms for E and f2 by the 
method we described above, we can write these in the desired form, namely 



E = P -e~ 2A + V -e~ iA log e~ 2r + P -e~ iA 
2 8 & 2 



i + V - + 1) - ^(1)) 



+ 



(3.64) 



r r z \ 6p z 



(3.65) 



Inserting these expansions into the expressions for the momenta, after taking into account 
the effect of the counterterm 



2 Vr 



2r 



(3.66) 



we obtain the two-point functions 



16 

(0(p)0(- P )) = ( 4^(1) - 4^(1 + «) + ^2 



<2?(P)0(-P)> 



\/3k 



2(0) b 



pV'(2ii(p)0(-p)> 



{Tij(p) T ki(-p)) 



TT 



n 



v3k 

ijklT, o 



/ = (0)V (Ti(p)T](-p))=0 



1 



(^(a + l)-V(l)) 



4. AdS-sliced Domain Walls 



AdS-sliced domain walls have also been studied in the literature p7L |2|, |2|, 111, P]. In this 
case the background is of the form 



ds 2 B = dr 2 + e 2A ^g ij (x)dx i dx i , $ = ^ B (r), 



(4.1) 
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where gij{x) is the metric of Euclidean AdSj with radius I and we have set the radius of 
the bulk AdS(i+i equal to l 6 . Inserting this ansatz into the bulk equations of motion leads 
to the following equations for A(r) and (pB( r ) 

A + dA 2 + ^-V^ B ) + i^ e ' 2A = 0, (4.3) 
a — 1 l z 

4> B + dA<p B -V'(0 B )=O. (4.4) 

Note that as I — * oo these reduce to the equations for fiat domain walls, as they should. 
From now on we set I 2 = 1. 

4.1 Janus solution 

A particularly interesting AdS-sliced domain wall solution is the dilaton domain wall so- 



lution of type IIB supergravity of [14] . This is a non-supersymmetric regular solution. 



When reduced to five dimensions, it solves the field equations of AdS gravity coupled to 



a massless scalar with a constant potential. Similar solutions exist in all dimensions [15]. 
These solutions are of particular interest because they enjoy non-perturbative stability 
for a broad class of deformations [15]. This strongly suggests that they should have a 
well-defined QFT dual. 

Setting V = — ■> the equation for the scalar field can be trivially integrated to 

give 

(j>B = ce~ dA , (4.5) 
where c is an arbitrary constant of integration. The remaining equations imply, 

A 2 = 1 - e~ 2A + be~ 2dA (4.6) 

2 2 

where b = g7gz/n • ^ ne geometry is non-singular provided the parameter b is within the 
range, 

1 ( d -i^ d - 1 



0< b < bo = ~ I — ) (4.7) 



One can obtain an implicit solution of (JO*) as 



\ A dA 

' o Vl - e~ 2A + bF 2 ^ 



6 In |15| a different convention was used: the radius of bulk AdSd+i and of the AdS^-slice were set equal 
to each other. 

7 Additional dilatonic deformations have been presented in |3C|]. 
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where Aq is the smallest zero of P(u) = bu d — u + 1, where u = e~ 2yi . This defines half 
of the geometry, i.e. the region with < r < oo. The other half is obtained by extending 
A(r) to negative values of r as an even function, A(—r) = A{r). 

We can obtain an explicit expression for the bulk metric by changing variables from r 
to u. Using 

A 2 = 1 - u + bu d 7 (4.9) 

we obtain 8 

du 2 1 

*b = 4u 2 (1 _ u + bud) + -%-(*)<fcW. (4.10) 

Note that if 6 = this is precisely the metric for AdS^+i in the ^cLS^-slicing parameter- 
ization. The range of the u-coordinate depends on the value of the parameter b, namely 
< u < u , where u Q > 1 with equality iff b = 0. We give the explicit form of 
function of b in appendix |A[ In this parameterization the two halves of the space, i.e. r > 
and r < 0, are not distinguished since u is an even function of r. In particular, the regions 
at r — ► ±oo are mapped to u = 0. 

We discuss in appendix ^ the conformal compactification of the solution. The confor- 
mal boundary consists of two half-spheres with angular excess joined along their equator 



1 14, 15]. We will refer to the joining equator as "corner". In order to calculate correlation 
functions of the dual field theory we need to write the Janus metric in the Fefferman- 
Graham (FG) form. Provided the boundary metric is smooth, this is always possible in a 
neighborhood of the boundary but the FG radial coordinate may in general not be valid 
far away from the boundary. In the present case, the boundary metric is smooth except 
for the presence of corners. We therefore except to be able to find a FG coordinates that 
are well defined in the neighborhood of the boundary except perhaps at the corner. 

In appendix [B] we construct the FG metric to all orders in b for the Janus geometry and 
determine the range of validity of the radial coordinate. We find that the FG coordinates 
are well defined everywhere in a neighborhood of the boundary except on the corner where 
the two half-spheres of the boundary meet. In particular, the FG metric takes the form 

ds 2 B = \ [dz 2 a + (1 + bc 3 (x) + 0{b 2 ))dz 2 d + (1 + 6c 4 (x) + 0{b 2 ))dzl] , (4.11) 

z o 

where x = Zd/z Q and z a , a = 1, . . . ,d — 1 are the standard transverse coordinates in the 
upper half plane parameterization of the AdS^ slice. The location of the corner is at 
Zd = 0. The functions c%(x) and c±(x) as well as the form of the FG metric to all orders 
in b are given in appendix [B[ As discussed there this coordinate system covers the region 



8 An equivalent form of this metric with A instead of u as a variable was found by C. Nunez (unpublished 
notes, July 2003). 
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\x\ > x = b/s/2 + 0(b 2 ), so z d = only when z Q = 0. In other words, this coordinate 
system does not cover a (radially extended) neighborhood of z d = 0. 
In this coordinate system the background scalar takes the form 

Mx) = o + cc 5 (x) + O(c 3 ), (4.12) 

where again c${x) is given in the appendix. It is significant to point out here that on the 
boundary, i.e. z Q = 0, the value of the scalar field is a step function in z^, namely 

<pB(z d ) =4> + sgn(z d )c, (4.13) 

which implies that the coupling of the dual operator is different on the two sides of the 
corner, or 'wall', at Zd = 0. These results are sufficient for calculating correlation functions, 
which we do in the next section. 

4.2 VEVs 

Now that we have determined the appropriate FG coordinate system we can carry out 
the algorithm developed in |J and evaluate the vevs of the stress tensor and the scalar 
operator dual to the dilaton, as well as, all two point functions using perturbation theory 
in c. The first step is to define the radial coordinate 9 r = — log z Q which is used as the 
'time' coordinate in the Hamiltonian formalism of ||. Due to the fact that the background 
depends also on the transverse space coordinates, a full asymptotic analysis is required to 
determine the covariant counterterms. We will not give these here but they are easily de- 
termined following the procedure in Q. Evaluating these counterterms on the background 
using the following expressions for the non-vanishing components of the Christoffel symbol 
and Ricci tensor: 

F d dd = \e%{x) + 0(b 2 ), r a bd = h -e^ 4 {x)5i + 0(b 2 ), (4.14) 

R dd = - { -^^e 2r c'i(x) + 0(b 2 ), R ad = R ab = 0(b 2 ), R= J±^4(x),+0(b 2 ) 

and adding them to the canonical momenta obtained directly by differentiating the back- 
ground fields w.r.t. r one obtains the following expressions for the vevs of the scalar 
operator and the stress tensor: 

(0)b = , (T}) B = 0. (4.15) 

Although the calculation has been done to leading order in c, it is not difficult to show 
that these results are in fact exact. The reason is that the coordinate transformation (|A.3[) 



9 This radial coordinate is different from the original radial coordinate in (4.8) but we hope this causes 
no confusion. 



- 19 - 



ensures that for every power of b there is a factor of z^ d ^ which means that higher order 
in b terms are subleading and do not survive when the regulator is removed. This can also 
be seen from the exact expressions for the Fefferman-Graham metric and scalar background 
given in appendix [B], which can be used to obtain the exact canonical momenta. Namely, 

K Bd = f 1 + 73^) 7 ' K% b =(u+ ^(l- u )(l-u + bu^ 5 a b (4.16) 

and 

cj) B {x) = sgn(x)cu d/ Vl - u. (4.17) 

One immediately sees that the vevs given above are in fact exact, as claimed. 

The form of the vacuum expectation values is the one required by the symmetries of 
the problem. As shown in [p3[ , the 1-point functions for a conformal field theory on a flat 
space with a boundary at Zd = (which breaks the conformal group from 0(l,d + 1) to 
0(l,d)) are precisely of the form ( 4.15|) . In the present case we consider the theory on 



both sides of the wall Zd = 0. The McAvity-Osborn result applies separately to the two 
regions, > and < 0, and it gives 

(°)B = A d ; *d>0, (°)b = 7%, *d<0. (4.18) 

\Zd\ a \zd\ d 

In the present case c\ = — C2 = c. 

These considerations suggest |0] that the dual field theory for d = 4 is Af = 4 SYM 
possibly coupled to non-super symmetric conformal matter localized at z^ = and with 
9ym being different on the two sides of the wall (similar suggestions can be formulated 
in all dimensions). This is consistent with the symmetries of the model: the presence of 
the defect breaks the symmetries to 0(1,4) (i.e. the (Euclidean) conformal group in three 
dimensions). It would be interesting to investigate whether there is a classical solution of 
N = 4 SYM coupled to such defect that can reproduce (|4.15| ), but we will not pursue this 
here. 10 

4.3 Two-point functions 

Since the leading correction to the AdSd+i metric is order c 2 , while the leading corrections 
to the (off diagonal) two-point functions are order c, we can take the background to be 
exactly AdS and consider linear fluctuations driven by a source Tj which is of order c. 

10 Note added: A precise proposal for the dual theory was recently made in J34j: they consider A/* = 4 
SYM theory on two half-spaces separated by a planar interface that contains no matter and with a different 
coupling constant coupled to specific operator closely related to the TV" = 4 Lagrangian density. The field 
theory computations in |g4j exactly agree (to the extend that they can be compared) with the holographic 
computations described in this and next subsection. 
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Decomposing the metric fluctuations as was done for flat domain walls above we derive the 
following equations for the irreducible components: 



d-l 



-D g e) = 2K 2 ny j T* (4.19) 

ij = 2 K 2 -±T kd+1 , (4.20) 
d k 

f = -2 K 2 —T kd+1 , (4.21) 

2n 2 f d+ld+1 - e~ 2r Df] , (4.22) 



: ,,o lj n S < LV (syhiPon ■ ihsypun ^djSd ) . (1.2:5) 

Only the first and the last equations need further analysis as the rest give immediately the 
momenta as functions of the linear sources. The responses for both the transverse traceless 
metric fluctuation and the scalar field fluctuation can be obtained using the massless scalar 
bulk-to-bulk propagator 

°d tdjpfj d+l_d 2 



G «> = 2^ F {"— •2 + h( I (424) 



which satisfies 

I 

7~9 



-U g G{£) = 5(z, w) = —5{z - w). (4.25) 



Here c d = T(d) / '(T(d/2)ir d / 2 ) and £ = 2z w /(z 2 + w 2 Q + {z- w) 2 ). As z Q ->■ 



G(0~^K d (w,z) (4.26) 



where 



K d (w, z) = c d ( 2 ™° ) (4.27) 

\w£ + (w - z) 1 ) 

is the well known bulk-to-boundary propagator. To complete the calculation then we need 
the source which is 

= d^Bd v cl) + dvfod^ - g^ s 4> B (j> + + W ( 4 - 28 ) 

Here we have used the fact that the background scalar is a function of x = z d /z Q which 
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implies z d d Zd (f)B = 4>b- With a little more algebra the source can be cast in the form 

f d+ld+1 =4> B U- e -^ d *J>) + °( 6 )> 

f jd+1 = <p B (djt - ^<W) + 0(b), 



Ti 



ij = 9B 



— (di d dj(f) + 8 jd di(t> - Sijd Zd (j)) - e r 4>5ij 
z d 



+ 0(b). (4.29) 



Using these sources and the above bulk-to-bulk propagator we can now evaluate the canon- 
ical momenta which give the one-point functions with linear sources. It is not difficult to 
show that no counterterms contribute to the order c terms of the momenta. It turns out to 
be easier to obtain the two-point functions from the canonical momentum of the graviton 
by differentiating w.r.t. the scalar source rather than from the scalar momentum and so we 
only consider the graviton momentum here. Of course both calculations should give iden- 
tical results and we show this explicitly in appendix [C], where we calculate the two-point 
functions from the scalar momentum. 

To obtain the the canonical momentum of the transverse traceless component of the 
graviton we note that the inhomogeneous solution to its equation of motion is 

4 = 2k 2 J d d+1 w^w)G{0^k l jfj:{w) (4.30) 

so that asymptotically 



e 



2K 2 e- dr / d^w^^Kdiw^^kjffiw). (4.31) 



Substituting the above expressions for the source into the canonical momenta and differ- 
entiating w.r.t. the scalar source we arrive at the two-point functions 

nY^T/^O^)) = -2(d + l)cU i d d jI(z, w), (4.32) 

^^(T*(z)0(w)) = ^ [{0(z)) B d k 5^{z- w)) , (4.33) 

K l k {Tf(z)0(w)) = ~ [(0(z)) B d k 5^(z- w)) , (4.34) 

(Tl(z)O(w))=0, (4.35) 
where the projection operators are acting on z and 

I(Z,W) = C\ J d d xdx X 2 d+1 — - x 2 )(d+3)/2 + (f _ md [x 2 + _ ^2]d " ( 436 ) 
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The last three correlators can be easily seen as a consequence of the Ward identities ( 2.15| ) 



and (2.15), which is a non-trivial consistency check of our calculation. In fact, since the 
vevs are exact, these two-point functions must also be exact in c, although our calculation 
of the two-point functions was done only to leading order in c. 

The computation of the remaining of the 2-point functions (Tij(x)Tki(y)} and (0(x)0(y)) 
requires an analysis to order c 2 . This computation is rather complex since the background 
metric receives a correction at this order. This means that we need to linearize the bulk 
field equations around the corrected solution ( |B.6| ). The latter, however, is inhomogeneous 
and this complicates the analysis. Nevertheless, the conformal invariance of the boundary 
theory completely determines the two point function of the scalar operator, while the two- 
point function of the stress tensor is determined up to a scalar function (except for d = 2 
where it is fully determined) It would be interesting to check that the holographic 

calculation reproduces these two-point functions as well. 

4.4 Janus two-point functions vs boundary CFT 

Let us now take a closer look at the structure of the Janus two-point functions. We would 
like to show that the 2-point functions are of the form required by conformal invariance for 
a CFT on a space with a wall at z d = 0. 11 The subgroup of the conformal group 0(1, d+1) 
that leaves z d = invariant is 0(l,d). This is precisely the isometry group of the Janus 
metric. McAvity and Osborn [|33| have given explicitly the form of this two-point function 
in such a CFT. It is given by 

(THz)0{w)) = -sgn(z d V - )g U) X% - -A , (4.37) 



(d_ 1)^/2 \&J \ 3 d »3^ 

where 12 

v 2 = -^—, i = -r—i s = z-w, (4.38) 
£ + 1 iz d w d 

and 



a 8 



Xi = z d -dii = v ( —Si -rii) , (4.39) 



2z_d 

with rii = Sid. The normalization of the two-point function is fixed by the normalization 



of the vev of the scalar operator [|33j| . It is clear from (4.18) that this normalization has 
opposite signs for z d > and z d < 0, which is the origin of the sgn(z d ) factor. This 
expression applies for z d w d > 0, i.e. both points on the same side of the wall, but not 



n We are grateful to the authors of for pointing us to the work of McAvity and Osborn and for 
prompting us to check that our calculation is consistent with their results. 

12 We use £ here to conform with the notation of This should not be confused with the argument of 
the bulk-to-bulk propagator used earlier. 
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for z d w d < 0. The holographic expression (|C.8[) applies to both cases, however. Under 
conformal transformations that leave the hyperplane z d = invariant we have 

-2 & Z d Wd Ar .s 

s -> n ^ w _ v z d ^——, w d ^——. 4.40 

It follows that £ is a conformal invariant while transforms as a vector. In particular, 
under inversion z — ► z/i^ 2 , w> — ► w/w 2 , 

Xi — > Iij(z)Xj, (4.41) 

where hj{z) = 6ij — 2^f. It is easy then to see that the two-point function given above 
transforms correctly under inversion, namely 

(T;(z')0(w')) = ^w^WliWiTftfOiw)), (4.42) 

One can show in general, using the fact that the background has the correct isometries, 



that the holographic 2-point functions transform as they should. Since the results of [33] 
follow from the same symmetries, this argument shows that our results are consistent with 
that of [B3]. It is, however, a rather non-trivial exercise to explicitly demonstrate that 



the correlator is of the form given in [33], mainly because of the integral representation of 
the transverse-traceless part of the correlator. The integral that appears in the transverse 
traceless part of the holographic two-point function is not easy to evaluate in general, and 
evaluating the projection operator acting on it is not straightforward either. This makes 
a direct comparison of the two results rather non-trivial. Instead, we will expand both 
results in a short distance expansion and compare them term by term. We do this for the 
first three orders in the expansion and we find complete agreement. 

To facilitate the comparison we first expand the above result of McAvity and Osborn. 
Of course, this expansion is valid only when z d w d > 0, which is also the condition for the 
validity of the McAvity-Osborn expression. After some algebra we get 

WW)) (d _ 1)7rd/2 \ 2wd \d+i \ g2 d °j 

1 / ■ ■ s*s 



. n l Sj + rinS 1 — ft ■ s5\ + (d — 2)n ■ s- 
2w d V 3 s 2 



1 

4 



~ (d(d - 2){n ■ s) 2 - {d + 2)s 2 ) S -^- + 2d(n ■ sjn^s^ + sVn, 

(4.43) 



(2w d ) 2 

l -{(d + 2)(n.s) 2 -s 2 )5^ + 0(s 3 )} 



The holographic result can be easily evaluated to this order too. First, using 

r(|) D i 

2(d-2)vr d / 2 fg2\ 



6 (d) ® = - ^/^L^ TT^f. ( 4 - 44 ) 
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we find that the longitudinal part of the holographic two-point function reproduces precisely 
the first two orders of the short distance expansion of the McAvity and Osborn result. The 
transverse traceless part is then evaluated by acting with the projection operator on 

and it reproduces exactly the third order term. Some details of this calculation are pre- 
sented in appendix Therefore, at least to this order in the short distance expansion, we 
have shown that the holographic two-point function is exactly what one expects for a CFT 
with a wall at Zd = 0. 



5. Conclusions 

We discuss in this paper the computation of correlation functions for holographic RG flows. 
The computation was done within the Hamiltonian framework developed in Q. A central 
point in our discussion is that the analysis is focused on the canonical momenta which are 
associated with regularized correlation functions and not the on-shell action. Furthermore, 
the renormalization procedure is set up such that one only computes the contribution 
of counterterms to correlators under discussion rather than first computing the general 
counterterms and then specializing. In particular, to renormalize n-point functions we 
only need to obtain asymptotic solutions to (n— l)-order in fluctuations. For instance, for 
2-point functions a linearized analysis is sufficient. This leads to significant reduction of 
labor compared to previous works. 

In the literature the analysis of fluctuations and of renormalization were often per- 
formed in different coordinate systems. This was due to the fact that renormalization 
was heavily based in the universal form of AdS asymptotics which required the use of a 
particular coordinate system, the Feffer man- Graham coordinate system. The fluctuation 
equation however is often more easily solvable in other coordinates. In order to combine the 
renormalization results with the solution of the fluctuation one needs the transformation 
between the two coordinate systems. This is straightforward to obtain, at least asymptot- 
ically, which is the only thing needed, but it adds to the complexity of the method. In 
this work we also avoid this complexity, since the asymptotic analysis is done directly at 
the level of the fluctuation equations. Moreover, it is not necessary to use the Fefferman- 
Graham coordinates anymore since the asymptotic expansion is now done by organizing 
the asymptotic solution in terms of eigenfunctions of the dilatation operator H. 

To illustrate the method we discussed both flat and AdS-sliced domain walls. We 
reduced the computation of 2-point functions that involve the stress energy tensor and 
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the operator dual to the active scalar for the most general flat domain walls driven by 
a "superpotential" W to the solution of two second order linear ODEs or (equivalently) 
to two first order non-linear ODEs. With the new method it is easy to carry out the 
renormalization in general (and we discussed how to do this) but it is even easier - almost 
trivial - to carry out the procedure in each specific case. For comparison purposes, we 
discussed all details for the two cases mostly studied in the literature, the GPPZ and 
CB flows. The new method is comparable in efficiency to the "old approach" where one 
extracts the correlator from the leading non-analytic part of the linearized solution, but it 
does not suffer from any of its drawbacks. 

As a new example, we discussed the computation of correlation functions of the Janus 
solution. This is a regular non-supersymmetric but stable AdS-sliced domain wall solution. 
The main difficulty in this example is that the boundary has a corner. However, we showed 
that there exists a Feffer man- Graham coordinate system which is well-defined everywhere 
in the neighborhood of the boundary except on the corner. This allows the vev's of the dual 
operators to be read off. We found that the expectation value of the stress-energy tensor 
is zero and the expectation value of the operator dual to the active scalar is non-zero. The 
expectation value is non-homogeneous and blows up at the location of the corner, which 
plays the role of a 'wall' in the boundary CFT. We computed then the 2-point functions 
that receive a contribution to leading order in c, i.e. {Tij(x)0(y)). Ward identities relate 
(part of) (Tij(x)O(y)} to the vev {(D). It is a nice check of both the value of the vev and 
of the 2-point function that the Ward identity is indeed satisfied. Both the vev and the 
2-point function we computed are of the form implied by conformal invariance for a CFT 
on R d with a wall at = 0. 

We discussed here only 2-point functions. Higher point functions have been discussed 
in the context of holographic renormalization in [J5], 35 1 . The new method can be straightfor- 



wardly applied in such cases, essentially trivializing the issue of renormalization. Another 
direction in which the new method holds promise in delivering new results is the case of 
solutions that are not asymptotically AdS such as the Klebanov-Strassler solution |3G[|. 
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Appendix 



A. Conformal compactification of the Janus solution 



In order to determine the conformal compactification of the Janus geometry we introduce 
a new radial coordinate 

z = A = ±\A -u + bu d . (A.l) 

This coordinate has the range — 1 < z < 1 for any value of b and the u = region is mapped 
to z = dbl. u can be determined as a function of z by solving the algebraic equation 



u-bu a 



1 



(A.2) 



as a power series in b. The relevant solution is the smallest real positive root which is given 
for arbitrary d and to all orders in b by 

r(nd+l) 



(d-l)-F(d-2) 



u(z;M) = (l- Z 2 )X; f 
d-l 



1 2 

d' d 



d 



n=0 

2 



(n + l)r(n(d - 1) + 2) 
3 d-2 d 



d d 

(d^T)^ 6(1 



(A.3) 



J2\<2-1 



d- 1' d- 1' ' d- 1' d- 1 
where P -F g is the generalized hypergeometric function. Note that the bound of the u 
coordinate mentioned above is just u a = u(0;b,d). It is not possible to express u(z) in 
terms of elementary functions except for the cases d = 2 and d = 3. We have respectively, 



u 



(*; 6,2) = i(l-Vl- 46(1-^) 



u(z; 6, 3) 



36 



sin 



— arcsin ( — V3b(l — z 2 ] 



(A.4) 



If we now write the (Euclidean) ^dS^-slice metric in global coordinates and set z 
sin#, the metric ( [1.10 ) becomes 



ds% 



u(sin 9) cos 2 A 



d\ 2 + cos 2 A 1 + (2d - l)6(cos 



,2, 



1-1 



+ 0(b 2 )j d6 2 + dr 2 + sin 2 AdO^_ 2 

(A.5) 



where < A < ir/2 and —tt/2 < 9 < ir/2. A few comments are in order here. First, the 
transformation ( |A.3 ) implies that every power of b in the coefficient of d9 2 comes with a 
factor of (cos 2 #) rf_1 and hence the metric inside the square brackets is non-singular for any 
9. Second, note that the (A, 9) part of the metric can be transformed into the standard 
metric on S 2 by introducing the angular coordinate 



li 



dz 



o \/u{z) (1 - bdu(z) d - r ) 



1 



1 3 



+ d - - bsin9F -,- -d;-;sin 2 +0(b z ). 



2' 2 



(A.6) 



This is precisely the angular coordinate introduced in [l5|] and it takes values in 
[— fioj/io], where \i > n/2 is given in equation (B.8) of [15|. Because of the excess an- 
gle the compact metric has a corner at A = tt/2, as is discussed in 
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B. Feffer man- Graham coordinates for Janus metric 



To construct the Fefferman-Graham metric we start with ( |4.10 ) and the coordinate trans- 
formation ( A. 3 ) and write the AdS d-s\ice metric in the upper-half plane coordinates. Then 

dz 2 



ds 2 B 



{l-z 2 ) 2 
1 



l + 2(d- 1)6(1 -z 2 ) d ~ v + 0(b 



+ 



l-z 2 



1 



1 - 6(1 - z 2 ) d ~ x + 0(6 2 ) -^(d~z 2 +dz 2 a ), 



where a = 1, . . . ,d — 1. For 6 = the coordinate transformation 

Zd - 



zl + z 2 



Z 2 o + 4 



(B.l) 



(B.2) 



brings this metric into the upper half plane metric with radial coordinate z . To determine 
the Fefferman-Graham form of the Janus metric we need to obtain appropriate 6-dependent 
corrections to this transformation. We can determine these as a Taylor series in 6 by 
introducing two arbitrary functions at each order in 6 and solving the differential equations 
that result by requiring that the transformed metric is of the Fefferman-Graham form. 
The unique transformation which ensures that the metric remains asymptotically AdS 
independent of 6 is to linear order in 6 

Zd 



% + 4 



+ bh(x) + 0(b 2 ), z = Jz 2 + z 2 + bz f 2 (x) + 0(6 



(B.3) 



where x = Zd/z Q and 
A 0*0 



X 



2 ( 1 + x 2)3/2 



1 'Yd) jr d F ( d ' d > d + l i-ia \ + 



2(d + l)x 2 ( d + 1 )' 



-F ld+l,d;d + 2; ^ + 



(1+x 2 )^ 1 



(B.4) 



f2(x) 



2V1 + 



1 + 2dx 2 ^ ( , , , 1 1 , 



1 



2(d+ l)x 2d 
The metric then takes the from 
1 



F[d+l,d;d + 2;-- j {1 + ^ 



1 



(B.5) 



ds% = — [dz 2 + (1 + 6c 3 (x) + 0{b 2 ))dz 2 d + (1 + 6c 4 (x) + 0{b 2 ))dz 2 a ] , (B.6) 



where 



^) = ^^nd,d-,d+i-,-±)-. 



x 2 ' 2{d + l)x 2 ^ d+1 ) 



' F(d+l,d;d + 2;-^) 



c ^ x) = ~^h 2dF{ - d,d]d+v, ~h 



X 2(l+ X l)d-V 

(B.7) 
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Note that the derivatives of these functions have a much simpler form: 



1 



x 3 (l + x 2 ) d 



-l-(2d-l)^ + 2 1 + ^ + 2 ^ 



21 



1+X 



2\2 



c' 4 (x) 



x(l + X 2 ) d ' 



(B. 



The metric in ( [B.6D is manifestly invariant under translations and rotations of the z a 
coordinates and scale transformations (the x coordinate is invariant under scale transfor- 
mations). The original metric (4-10) however was invariant under the larger group 0(l,d) 
associated with the AdS slice metric. We now show that the metric ( |B.6| ) is also invariant 
under a discrete inversion isometry to order b which enhances the isometry group to the 
full 0(l,d). Actually, we will see that the inversion symmetry can be used to obtain the 
Feffer man- Graham form of the metric to all orders in b. 

Let us write the AdS slice metric in ( 4.10| ) in the upper half plane coordinates so that 

du 2 



ds 2 B 



Au 2 {\ — u + bu d ) uz 2 
This form is invariant under the discrete isometry 



+ -=s{dxl + dz a dz a ). 



(B.9) 



~2 i 2 : 
4 + Za 



zl + z 2 ' 



(B.10) 



We now bring this metric into the Fefferman- Graham form by means of a coordinate 
transformation 13 

z = s(x;b), z = z Q t(x;b), (B.H) 

where x = Zd/z a . We point out that this is precisely the form of the coordinate transforma- 
tion ( |B.3| ), but we now treat the 6-dependence non-perturbatively. This allows us to express 
the above discrete isometry in terms of the Fefferman-Graham coordinates z^ = (z Q , z a , z&). 
We find 

(B.12) 



z 2 t(x; b) 2 + z 2 



Now, the Fefferman-Graham metric (B.6) takes the form 



ds 



B 



[dz 2 + X(x; V)dz\ + fj,(x; b)dz a ] 



(B.13) 



The requirement that this is invariant under inversion uniquely determines the functions 
\(x; b) and fj,(x;b) in terms of t(x; b). Namely we find the exact FG metric 



ds 2 -1 
as B - 



dzl + 



d x t 



x(t — xd x t) 



dzl + 



t(t — xd x t) 



dzi 



(B.14) 



Note that u = u(z; b) is given in eq. (A. 3) 
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Requiring further that this is equal to the Janus metric above uniquely fixes the trans- 
formation functions s(x;b) and t{x;b). In particular we obtain the system of coupled 
equations 

(B.15) 



u = i _ x ^l } (d x s) 2 = —u 2 (l - u)(l - bdu d ~ 1 ) 2 , 



dxt 
t ' 

where u = u{s{x)) is given by ( |A.3| ) . 

In order to solve these equations we use ( |A.1| ) to trade s(x) for u(x) in the second 
equation, which gives 



(B.16) 



The sign and the integration constant are chosen so that u(x) ~ 1/x 2 as x — > oo, indepen- 
dent of b. Unfortunately it seems rather difficult to do this integral explicitly for arbitrary 
dimension d. Instead, one can expand the integrand in b and integrate term by term. This 
gives 




1 + x 2 2(2(1 +x 2 ) rf + 2 ^ ' 1 + x 2 ' 

The transformation functions s(x;6) and t(x;b) are now determined from 



(B.17) 



s(x) 


= 1 - u(x) +6ii(x) d 


i(x) 




"l du' / 1-n' \ 1/2 ~ 




= exp 


2Ju(x) u> [l-u' + bu'dj 





(B.18) 



Inserting the above expansion for u{x) we reproduce (after some manipulation of the hy- 



pergeometric functions) precisely the coordinate transformation (B.3) 
Moreover, inserting ( B.lTj ) in 




gives 

where 4> Q is a constant and 

C5O) 



<j> B (x) = (j) + cc 5 (x) + C(c 3 ), 



x /l d 3 x 2 
^ V2 ! ~ 2' 2' 1 + x 2 



vT+x 



(B.19) 



(B.20) 



(B.21) 
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Again, this has a simple derivative: 

^ = {1 + J)(W (B.22) 

Notice that as z Q — ► with all other coordinates fixed (i.e. as we approach the conformal 
boundary) c§(x) = sgn(z^) while higher order terms do not contribute. So at the boundary 

4>B{zd) = 4>o + sgn(z d )c. (B.23) 
This implies that the coupling of the dual operator is different on the two sides of the wall. 



Finally, let us examine the range of validity of the coordinate transformation ( B.ll ) . 
The Jacobian of the transformation is equal to J = td x s. Now J = implies d x s = 
since t(x) is positive definite, as can be seen from ( |B.18| ). It follows that the coordinate 
transformation breaks down at it = 1. Note that the zero of (1 — bdu' 1 ' 1 ) occurs at 
u = 1/ (bd) l ^ d ~ 1 ^ > 1, where the inequality follows from ( |4.7| ). We conclude that the 
Feffer man- Graham coordinates are valid in the range < u < 1 although, in general 
< u < u Q with u Q > 1. Recall that in general the Fefferman-Graham coordinate system 
is only guaranteed to exist in a neighborhood of the boundary, and here we see an explicit 
illustration of this. 



Recall that the Fefferman-Graham coordinate system [31] is obtained as follows (see 
section 3 of [|j for a review). One considers Gaussian normal coordinates centered at the 
boundary and the radial coordinate is identified with the affine parameter of the geodesies 
emanating perpendicularly from the boundary. Clearly the region of validity of this coordi- 
nate system depends on the behavior of the radial geodesies. We therefore need to analyze 
such geodesies, and we will do this in the (u, z Q , z a ) coordinate system which is well-defined 
everywhere. 

One easily shows that there are geodesies with z a constant. The geodesic equations 
for the remaining coordinates lead to the following two equations 

d log z„ 

= a lU , (B.24) 



+ — — J -)u 2 + 2a 2 1 u 2 (l-u + bu d ) = 0, (B.25) 



u 2(1 -u + bu d , / 

where a\ is an integration constant. If a\ ^ the second equation can be integrated once 
to get 

it = ±2a 1 u^ (a 2 -u)(l-u + bu d ) (B.26) 
for some constant < a 2 < u a . If a\ = one gets instead 



u 



±a 3 ii\/l - u + bu d , (B.27) 
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Figure 1: The qualitatively different radial geodesies of the Janus geometry: (i) a 2 < u Q with 
a\ > (solid arrow) or a% < (broken arrow), (ii) 02 = u with a\ > (solid arrow) or a\ < 
(broken arrow) and (iii) a\ = 0. These qualitative features are insensitive to the value of b. 



where 03 is again a constant. Now, depending on the values of the parameters a± and 02, 
we can identify three qualitatively different types of geodesies as shown in fig.|l[ 

Consider now the radial geodesies defined by = z a = in the Fefferman-Graham 
coordinates, where the dot stands for the derivative w.r.t. the afHne parameter r = \ogz Q . 
Since z^ = constant along these geodesies we will take r = \og(z /\zd\) = — log \ x\ for later 



convenience. The transformation (B.ll) immediately gives 

dlogz 



dr 



u, (B.28) 
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Figure 2: A radial geodesic in the Fefferman-Graham coordinates (left) is not denned beyond 
u = 1 due to the failure of the coordinate system at this point. In the coordinate system which 
extends beyond u = 1 this geodesic bounces back at u — 1. Only the branch before the bounce 
corresponds to the geodesic on the left. 



while (B.16) implies 



u = 2u^(l-u)(l-u + bu d ). (B.29) 

The Fefferman-Graham radial geodesies therefore correspond to radial geodesies with a\ = 
a,2 = 1. In particular, they are geodesies of type (i) if b > but they are type (ii) if 
6 = 0. This is an important qualitative difference between the FG coordinates for the 
Janus geometry and pure AdS. This is in fact why the FG coordinates cover the whole of 
AdS but only part of the Janus geometry. 

It is now clear why the FG coordinate system for b > breaks down at u = 1. Namely, 
the radial FG coordinate corresponds to geodesies which do not reach beyond u = 1. If one 
continues to affine parameter values greater than r*, where u(r*) = 1, the geodesies bounce 
back and they cannot be used to define a coordinate system since they doubly cover the 
region u < 1 as is shown in fig.Q. Therefore the FG coordinates are well-defined for affine 
parameter values r < t*. This means that \x\ = e~ T must be bounded below. Another way 
to see this is to observe that ( |B.16| ) implies that x 2 is a monotonically decreasing function 
of u. Hence the upper bound u < 1 on u implies a lower bound on x 2 . Setting u = 1 in 
(B.17) and solving for x to leading order in b we find 14 \x\ > x Q = b/V2 + 0(b 2 ). Therefore 



x = is not part of the manifold and hence the metric ( |B.6| ) (and ( |B . 1 4| ) ) is non- singular in 
the region it is well-defined. (To cover the entire spacetime one would have to use another 



14 One must be careful since the hypergeometric function F(d,2; d + 1; ^ ) is singular at x 2 = 0. See 
eq. 15.3.12 in pi. 
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coordinate patch that covers the deep interior region 1 < u < u Q , but this is irrelevant 
for our holographic computations.) Notice that the bound on x translates into \z<i\ > x z 
and so only at the boundary z Q — ► does Zd cover the entire real line. More crucially, 
the bound z Q < \zd\/x Q means that the radial coordinate, z D , is well-defined everywhere 
except at z^ = 0, which precisely corresponds to the corner where the two-halves of the 
boundary meet. This can also be seen directly from the properties of type (i) geodesies. 
Since z Q — * \za\ as r — > —oo with z^ constant along the geodesies, we have 

z (t) = \zd\ exp / (1t'u(t'). (B.30) 

J — oo 

So for b > the FG geodesies hit again the boundary u = at z Q = \za\ct, where a = 
exp L , du At Zd = these geodesies become degenerate and they stay along 

KJ0 yj(l-u){l-u+bu d ) as, 6 

z Q = 0. Since the entire z a = subspace is mapped to (z a , Z&) = (0, 0) in the FG coordinates, 
the FG geodesies do not leave the origin once z^ = 0. Hence, the FG radial coordinate is 
not defined at z^ = 0. 

Finally let us discuss the possibility to use the geodesies of type (ii) in order to define FG 
coordinates. This is also a natural choice as these geodesies are the obvious generalization 
of the pure AdS case. Following radial type (ii) geodesies for the Janus geometry, the exact 
FG metric is 

.2 



= d 4 + U ° {U °: U) dzj + 4 (B.31) 
zi uzj uzi 



where 



and 



log (^vUf - dU (B.32) 

[4 J u^{u -u){l-u + bud) 



zl = z 2 d u ° exp / = (B.33) 

o \J{u — u)(l — u + bu d ) 



Asymptotically, i.e. as z Q — > 0, 

1 



~ ry 



dz 2 + (z 2 d )^o \dz 2 + dz 2 a ) 



(B.34) 



where Zd = z d °. So when u Q = 1, i.e. the pure AdS case, we get the standard result but for 
the Janus solution these geodesies lead to a non-flat representative of the conformally flat 
conformal structure. One can perform the additional change of variables given in (8)- (9) 
of |j?J to change the representative to the flat metric. Notice however that the conformal 
factor is singular at Zd = so the corresponding coordinate transformation is singular 
there. We thus arrive again at the conclusion that the FG coordinates are not well defined 
at the corner. 
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We have therefore determined the exact form of the Fefferman-Graham metric for the 
Janus geometry and we have shown that is well-defined everywhere except on the defect 
where the two half-boundaries are joined and it is non-singular where it is defined. More- 
over, we have shown that the FG metric possesses an inversion isometry which enhances 
the isometry group to the full 0(l,d) isometry group of the original Janus metric. This 
is reflected in the fact that the holographic calculation gives a zero vev for the stress ten- 
sor, which is consistent with a boundary QFT invariant under conformal transformations 
leaving the plane z d = invariant. 



C. Two-point functions for Janus from scalar momentum 

Here we give some details of the calculation of the two-point functions for the Janus back- 
ground based on the scalar equation of motion ( 4.23[ ). The non-trivial part of the calcula- 
tion consists in casting the source in a form which significantly reduces the amount of work 
required. To this end we again use the fact that 4>b{x) is a function of x = Zd/z Q only and 
it satisfies 

4>B + d<j) B + e- 2r U(t)B = 0(b) 



to write 



1 



dj(p B (x) = 5 jd —<j)B, didjcpB = -(d+ l)5 id 5 jd 



zi + zt 



+ 0(b). 



(C.l) 



(C.2) 



Decomposing S*j into irreducible components then ( 4.23j ) becomes 

1 



-a 



-cz a di 



(1 + X 2 ) 2 



+- 



cx 



(1 + x< 



d+l 



2(d-l 



-□/ + 



(d + l) 
1 + x 2 



2d ZH e d + 



d 



d-l\d □ 



d 2 



□ 



Z d 



De d -d z J + -d Zd S 



+ 0(b). 

Quite remarkably, this can be cast in the form 



where 



-a 



cz" 



-cz Q di 



(l + x 2 



d+l 



n g J^ + 0(b), 



1 d 7 

a + e d + ^S + 



1 



J^Zd \ r 

2D" ' 2(d-l) \ Zd ~ TT 1 1 



and a is a constant. Hence, the inhomogeneous solution is 

j„-c i d d+i w^^)G(m 



(1 + x s 



d+ 



T e d 



(C.3) 



(C.4) 



(C.5) 



(C.6) 
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where e d is given by the zero-order solution 

e){z) = f d d yK d (z,y)e {0) )(y). (C.7) 

This expression for (f> immediately gives the canonical momentum from which we obtain 
the two-point function by differentiating w.r.t. S'^j- 15 The result is 



{TAz)0{w)) = -2(d + l)cn i / i /(z,uJ) 



cd 



\w d \ d+1 



5 {d \z-w). 
(C.8) 



It is a straightforward exercise to verify that this is equivalent to the two-point functions 
given above, as calculated from the graviton momentum. 

D. Short distance expansion of the holographic two-point function (Ti(z)0(w)) 

For the convenience of the reader we will give here the essential steps required to evaluate 
the short distance expansion of the transverse traceless part of the holographic two-point 
function {T l -{z)0{w)), namely W d d jl(z, w), where I(z,w) is given in (|4.36j ). 16 

First, after a shift and rescaling of the integration variables, I(z, w) can be written as 

1^)=-^— r dXoX ^ [d d x wd+ ^ xd 1 1 



^) d/2 - 1 JO ° J ft? + ( Wd + \8\ Xd )^ K + ^ VI + (X - 8)*f 

(D.l) 

where s = s/s 1 . This form is suitable for a short distance expansion in | s| . Each term in the 
expansion can be explicitly evaluated using the standard Feynman parameters technique. 
The result to leading order is given in ( 4.45| ). 



To evaluate the projection operator on this expression we use the fact that 

! ' ;n7^-T> (D.2) 



(a 2 )* 2(a-l)(d-2a) [s 2 )^ 1 ' 



ir, 



Note that 



1 (O(x'))+5^(x,x')g {0)ij (x)(O{x')}, 



so if one starts the computation from the scalar 1-point functions, one should remember to include the 
contact term given above to obtain the full expression. 

16 Incidentally, this integral transforms under inversion as I(z',w') = z 2d w 2d I{z,w) and hence it must be 



of the form f{v)/{s 2 ) i for some function f(v), where v is defined in ( [f.38| ). However, we have not succeeded 
in determining this function so far. 
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for any power a ^ d/2 in order to cancel the 1/D factors in the projection operator. It is 
then straightforward to evaluate the derivatives in the numerator of the projection operator 
to obtain the short distance expansion of the transverse traceless part. The result is given 
in section 4.4. 
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